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Abstract 

We present two classes of nonequilibrium models with critical behavior. Each model is 
characterized by an integer q > 1, and is defined on configurations of q- valued spins on 
regular lattices. The definitions of the models are very similar to the updating rules in 
Wolff's algorithm for the Potts model, but both classes break detailed balance, except 
for q = 2 and q = oo. In the first case both models reduce to the Ising model, while 
one of them reduces to percolation (more precisely, to the general epidemic process) 
for q = oo. Locations of the critical point and critical exponents are estimated in 2 
dimensions. 

PACS numbers: 05.20, 05.40, 05.50, 36.20.E, 61.40, 64.70 
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It is by now well understood that nonequilibrium models can show critical behavior 
very similar to second order phase transitions in equilibrium systems. Indeed, 
dropping the requirement of detailed balance gives additional freedom, and one 
can observe a much richer spectrum of possibilities. These include in particular 
critical behavior for generic values of the control parameters, i.e. systems which 
do not require the control parameter to be set to any non-natural value to be 
critical. Examples for this are diffusion limited aggregation [p], |2| and various 
models of self-organized criticality f|, |5| . 

A problem which is more common in nonequilibrium models than in equilibrium 
cases is that it is often not clear how large are the universality classes. Thus, 
e.g., it is not clear yet whether different versions of the sand pile model || || 
are in the same universality class or not. An other much discussed example is 
provided by models with an absorbing state. The prototype of such models is 
directed percolation, interpreted as an epidemic process without immunization. 
By now it seems clear that a model for heterogeneous catalysis proposed some 
years ago by Ziff et al J7| does belong to the same universality class, in spite of 
original numerical indications of the contrary ||. But there exist still a number 
of similar models (some with degenerate absorbing states) for which conflicting 
evidence has been reported as concerns universality. 

Another problem which is more frequent in nonequilibrium cases is that there 
exists no good field theoretic treatment. The best known example for this is 
diffusion limited aggregation which is still not amenable to such a treatment. This 



does not mean that one cannot write down "hamiltonians" || [L0| or path integrals 
IT ], but rather that no good perturbation expansions are available which then 
could be resummed by renormalization group methods. Other models for which 
this applies are the sand pile model [O] and "true" self avoiding walks |T3], for 



which the field theories are intrinsically non-renormalizable [jig] . 

In the present note we want to present two new classes of nonequilibrium critical 
phenomena. They are superficially similar to the Potts model fl6| , but they are 
not defined in terms of hamiltonians. Instead, they are defined algorithmically, 
by stating the rules for updating the state of a system. This is similar e.g. to the 
sand pile model. But in contrast to that, our models show "conventional" critical 
behavior in the sense that they involve control parameters, the critical behavior 
being observed only for special values. On the other hand, we have not been able 
yet to apply field theoretic methods to them. 

Both class of models involve "spins" with q possible orientations ("colors"), and 
are symmetric under cyclic permutations of the colors. In addition, the first model 
is invariant under the exchange of any two colors, just as the Potts model is. The 
rules for updating a configuration are very similar to Wolff's single cluster variant 
TTfl of the Swendsen-Wang dynamics [|18[] of the Potts model. Let us thus recall 



Wolff's algorithm as applied to the g-state Pots model. 
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Assume we have a configuration of spins on a hypercubic d- dimensional lattice, 
with spin values Sj, i G Z d . The evolution is defined by sequential updatings of 
randomly selected clusters. To make one update, we proceed as follows: 

(1) We pick a random color s G [0, 1, . . . q — 1]. 

(2) Pick a random site i G Z d . If Sj = s, we do not do anything and proceed to 
the next update. 

(3) Otherwise, we build a bond percolation cluster on the subset of sites j which 
have Sj = s, and which are connected to i. 

(4) After this is done, we flip this cluster, i.e. we change all its Sj to s, and proceed 
to the next update. 

The control parameter in this model is the probability p for connecting sites in 
the percolation process in step (3). It is related to the interaction strength K in 
the Potts model, e~^ H = exp(if J2 {itj) (5 SiiSj - 1)), by p = 1 - e~ K @. 

Our first class of models ("class A") is obtained from this algorithm by keeping 
steps (1), (2) and (4), but replacing step (3) by 

(3A) Build a bond percolation cluster on the sites j which are connected with site 
i and which have Sj ^ s. 

Again this is done with probability p for each bond. But for this model we were 
not able to prove detailed balance for general values of q, and we thus cannot 
relate p with a coupling strength in any hamiltonian. 

In the second class of models ("class B") we skip step (1), keep step (3), and 
replace the other steps by 

(2B) Pick a random site i G Z d . 

(4B) Flip all spins of the percolation cluster into color Sj + 1 modg, i.e. rotate 
the entire cluster in color space by an angle 27r/g. 

In this model the braking of detailed balance is evident. 

The simplest case is q = 2. In this case, the Potts model goes over into the Ising 
model if we perform the trivial relabeling s = 0— >s = — 1. It is easy to see that 
in this case also models A and B are equivalent to the Ising model. 

Another simple limit is q = oo. In this limit, Sj will be different from s with 
probability 1, and thus the cluster built in step (3A) is just an ordinary bond per- 
colation cluster on the entire lattice. Thus model A is equivalent to the spreading 
of bond percolation according to the general epidemic process |19|] or the Leath 
algorithm . 

The q — > oo limit of model B is rather different. It is easy to see that there is no 
ordering for any finite p in this limit. For finite q and p, there is a balance between 
ordering due to the fact that two sites i, k which had originally Sj = s^ — 1 might get 
the same colors by flipping j, and disordering because a coherent cluster is broken 
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up. For q — > oo, the ordering disappears, and thus the transition point has to shift 
to p c — > 1. For q = oo all clusters are trivial (just one site) for all p, suggesting 
that the transition turns into first order (all scaling law amplitudes — > 0) for 
q — > oo. Notice that this conclusion rests on the fact that the fractal dimension 
of percolation clusters at p < p c is zero, whence the ordering by each cluster 
flip effects only a vanishing fraction of sites, and a random initial configuration 
remains essentially random. This is not so on finite lattices (in particular at 
d — 2), since there the codimension of percolation clusters is very small, and 
hence typical cluster flips involve large parts of lattice. Thus we expect very large 
finite size effects in model B even for finite values of q. 

Let us next discuss mean field theory. In both models, this implies that a spin not 
yet checked during the build-up of a cluster has the same probability 1/q to have 
any of the q possible colors. In a strict mean field treatment one would assume 
the same also for sites which had already been visited during the present cluster 
evolution. We shall in the following discuss a more realistic "hybrid" mean field 
ansatz where we do not make the latter approximation. In this case both models 
reduce to mixed site-bond percolation. In model A we have bonds established 
with probability p, and sites with probability 1 — 1/q. In model B, the bond 
probability is the same, but the site probability is 1/q. This implies correctly that 
model A turns into bond percolation for q — > oo. But it predicts a qualitatively 
wrong behavior for model B, as it would imply that model B does not become 
critical for q > l/p c , S ite, where p c , S ite is the threshold for site percolation on the 
same lattice. This is obviously a weakness of our mean field ansatz, as it would 
mean that even the Ising model (q = 2) is non-critical in 2 dimensions. The same 
problem would occur for the strict mean field treatment mentioned above, and for 
the original Swensen- Wang- Wolff model with q > 2. 

In following we shall report results from simulations on 2-d square lattices of size 
N x N with periodic boundary conditions. We used both depth-first [El], [22] and 



width-first ("Leath" |20[) algorithms for building the percolation clusters. While 
the first are somewhat simpler when implemented by recursive function calls, the 
temporal behavior is more natural in the latter as it corresponds to epidemic like 
spreading [l||. The critical point and the distribution of cluster sizes (and of 
cluster radii) is identical with both algorithms. In all sufficient number of 

cluster flips during the initial transient stage was discarded, i.e. all data reported 
below refer to the stationary state. 

In order to obtain p c for model class A, we used two different procedures. In 
the first (finite-size scaling) we used relatively small lattices (up to N = 1024), 
and determined p c ^ from the requirement that the cluster size distribution Pn(s) 
shows the longest scaling region for this value of p. In fig.l this is illustrated for 
q = 3 and iV = 1024. Values of p c ^ resulting from this and analogous plots are 
plotted against 1/N in fig.2. We see a straight line for N > 200 which extrapolates 
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to 

Vc = 0.5330 ± 0.0004 (1) 

In the second procedure, we used much larger lattices (up to 8192 x 8192) and 
sufficiently small values of p so that all cluster diameters were << N and finite- 
size effects were negligible. In fig. 3, resulting average cluster sizes (s) are plotted 
against p — p for three different values of p . Since we expect 

(s}~(Pc-p)-". (2) 

we expect a straight line at p = p c . This is indeed observed, with the same p c as 
above and with 7 = 1.65 ± 0.05. The error is mainly systematic due to a slight 
but significant curvature in fig. 3. 

Similar deviations form scaling behavior are seen in the average cluster evolution 
times and in the size distribution P(s). The former satisfies 

(T)~( Pc -p)- 5 , 5^0.85 (3) 

(see fig. 4), but there are rather strong deviations from a pure scaling law. 

From fig. 1 we see that the distributions Pn(s) on finite lattices at p = p c ^ fulfill 
roughly a scaling law 

P N (s)~l/s (4) 
for s << N 2 , i.e. in the usual notation [E31 we have r ~ 2. Results for larger 



lattices and p < p c are shown in fig. 5. We see there again roughly P(s) ~ 1/s for 
s < (p c — p)^ 1 (notice that bin sizes are oc s in fig. 5, in contrast to fixed bin sizes 
in fig.l). But this is superimposed by oscillations whose amplitude increases as 
we approach the critical point. The latter is very different from the Ising (q = 2) 
and percolation (q = 00) limits, but was also observed for q = 4 and for the time 
distribution for q = 3 and q — 4. On the other hand, model B leads to cluster size 
distributions which are qualitatively as in percolation and in the Ising model, see 
fig.6. 

Results for p c and for the critical exponent 7 are summarized in table 1, for both 
model classes. For class A, we see that p c interpolates monotonically between 
the Ising and percolation limits. The latter is approached very quickly. On the 
other hand, 7 seems first to fall below the Ising value 7/4, and rises only very 
slowly towards the percolation limit 43/18 = 2.388 for q — ► 00. For class B we 
see that there is a phase transition for all values of q with p c (q) — > 1 very slowly 
for q —>■ 00. The amplitudes in the scaling laws for class B seem to tend towards 
zero for q — > 00 (the clusters become very small except when p is very near p c ), 
in agreement with the arguments given above. We do not quote exponents 5 
because they are very strongly affected by corrections to scaling, and we do not 
quote exponents r for class B for the same reason. 

Superficially, our models resemble the Potts model. There, one has a first order 
transition for q > 4. One might thus wonder whether there is some range in q in 
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which one of our models has also a first order transition. We have not found any 
indication for that. 



In summary, we have in this paper presented evidence for two new classes of non- 
equilibrium critical phenomena. Both classes use Potts spins, and are defined via 
cluster evolution rules similar to those of the Wolff dynamics for the Potts model. 
The most interesting aspect of the first class is that it interpolates between the 
Ising model (q = 2) and percolation (q = oo). For this class we have not been 
able to find any observable which shows broken detailed balance. We thus cannot 
exclude that there exists for this model also an equilibrium version, possibly with 
long range interactions (otherwise it would be hard to understand why these 
models have not been found before). 

In class B, breaking of detailed balance is evident. Indeed, the transition in class 
B is from desynchronized cyclic behavior to synchronized one by formation of a 
percolating cluster of sites with common phase. Synchronization in spatially ex- 
tended models was studied recently in [^J and P5[ . In |24| it was argued that no 



synchronization can occur in stochastic discrete models with > 2 states. Contin- 



uous models were studied in [2S| where a connection was established with surface 
roughening in Kardar-Parisi- Zhang growth. The latter suggests that no phase 
transition with long range order occurs in d = 2 for such models either. These 
predictions are in clear contrast to our findings. This reflects basic differences 
in the dynamics underlying model B and the systems studied in p4| , p5| . While 
it is assumed in the latter that the phase proceeds rather uniformly with locally 
independent fluctuations, the phase progress in the ordered phase of model B is 
dominated by rigid flippings of entire large clusters. 

Finally, we have neither done any simulations in d > 2, nor have we attempted a 
field theoretic treatment. Using methods of |9|, 10], [TT[] it should not be difficult 



to set up the field theory. But it might be much harder to find a systematic 
perturbation expansion. 

Part of this work was done when one of us (P.G.) was visiting Rome University. 
He wants to thank A. Vulpiani and G. Paladin for their generous hospitality and 
for numerous enlightening discussions. 
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Table 1: 





= 2 




o = 4 


a = 8 


q = qo 


model A: 












Pc 


0.5858 


0.5330 


0.5170 


0.5029 


0.5 


7 


1.75 


1.65 


1.79 


1.88 


2.388 


T 


2 


~ 9 


~ 9 




2.055 


model B: 












Pc 


0.5858 


0.633 


0.6667 


0.7393 


1 


1 


1.75 


1.35 


1.21 


0.71 
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Figure Caption: 

Fig.l: Cluster size distributions for q = 3, lattice size N = 1024, and 3 different 
values of p: p = 0.53175(0), p = 0.532 (boxes), and p = 0.53225(A). Data were 
binned into bins of size As = 10 4 . The straight line has slope —0.9. 

Fig. 2: Finite size critical values p C: N versus l/N for q = 3. For large N the data 
seem to follow a straight line, the extrapolation of which to l/N = gives p c 
(dashed line). 

Fig. 3: Log-log plot of average cluster size for q = 3 against p — p , for p = 
0.5325, 0.533 and 0.5335. A roughly straight line is seen for po = p c ps 0.533, but 
notice that there are small but significant deviations. 

Fig. 4: Same as fig. 3, but for the average cluster life time. 

Fig. 5: Cluster size distributions (q = 3) for subcritical p on practically infinite 
lattices. Values of p are 0.5292, 0.531, and 0.5318. Data are binned into bins 
[s,2s-l]. 

Fig. 6: Similar as fig. 5, but for model B with q = 3. Values of p are 0.61, 0.63, 
and 0.6326. 
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